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Abstract. - Upper critical field, Hc2, in quasi-lD superconductors is investigated by the weak 
coupling renormalization group technique. It is shown that Hc2 greatly exceeds not only the 
Pauli limit, but also the conventional paramagnetic limit of the Flude-Ferrell-Larkin-Ovchinnikov 
(FFLO) state. This increase is mainly due to quasi-lD fluctuations effect as triggered by interfer- 
ence between unconventional superconductivity and density-wave instabilities. Our results give a 
novel viewpoint on the large Hc2 observed in TMTSF-salts in terms of a d-wave FFLO state that 
is predicted to be verified by the Hc2 measurements under pressure. 
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Introduction. — Thirty years after the discovery of 
Organic superconductivity in the quasi-one-dimensional 
(quasi-lD) molecular conductor (TMTSF)2PF6 [1], the 
fiature of this phase has not been entirely clarified in spite 
of intensive experimental and theoretical efforts devoted 
to this end. An enduring problem concerns the properties 
of superconductivity under magnetic field, in particular 
the issue of the upper critical field Hc2 in these materi- 
als, which is the topic of the present letter. The value of 
Hc2 determined for instance by resistivity measurements, 
drastically exceeds the so-called Clogston or Pauli limit, 
Hp. The excess occurs when the magnetic field is oriented 
parallel to the plane of (TMTSF)2J^ {X =PYq [2, CIO4 
P]), that is when pair breaking orbital effect is virtually 
quenched by anisotropy. There are basically two possibil- 
ities that have been put forward to explain the excess of 
Hc2'- the existence of a spin-triplet pairing or the presence 
of a Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. 

The constant spin susceptibility, Xj as obtained from 
the first Knight shift measurements across for X =PF6 
under pressure gave support to the former scenario for 
triplet superconductivity [3]. On the other hand, more re- 
cent NMR results on the ambient pressure superconductor 
X =C104, display a decreasing x below Tc, at least in low 
fields, which is compatible with singlet rather than triplet 



pairing [5]. As for the FFLO state, an estimation of the 
paramagnetic limit, 11^2^^, puts it lower than the ifc2 ob- 
served in (TMTSF)2X [S], which would make this scenario 
unlikely on a theoretical basis. However, recent detailed 
resistivity measurements on X =C104 compound using 
rotating field, display the signature of non-uniform super- 
conductivity at high fields, suggesting that a FFLO state 
could be realized [315]. The amplitude of Hc2 determined 
by the bulk measurements such as the nuclear relaxation 
rate [5], magneto-torque [S] and specific heat [TU], is lower 
than that of resistivity measurements, indicating the exis- 
tence of a high field interval H^f^ < H < i/^'f'", which 
we term as a "transient region" , where incomplete super- 
conductivity is realized. 

On theoretical side, the possibility of triplet supercon- 
ductivity in an interacting quasi- ID electron system has 
been investigated on a microscopic basis in various situa- 
tions. In zero field, /-wave triplet pairing has been shown 
to compete with d-wave singlet pairing, when long-range 
intrachain interactions become sufficiently large [TltiM] . 
The competition was also shown to emerge when inter- 
chain repulsive interactions are finite, albeit small |15[|16j : 
or as the amplitude of on-site interaction is huge |17| . 
A field-induced singlet to triplet pairing crossover has 
been also proposed to occur on phenomenological grounds 
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[T^14lll8j . A similar transition under field has been shown 
to take place microscopically using mean-field theory |19] , 
or when in a similar framework, long-range interactions 
are included [20H22] . The H-T phase diagram thus ob- 
tained, however, departs from observation obtained for the 
bulk [TllHllin]; it would also give a temperature indepen- 
dent spin susceptibility in the metallic state, at variance 
with experiments ■ 

The FFLO state in the quasi-lD geometry has been 
examined in great details in the framework of the BCS 
mean field theory, which predicts its existence in the region 
of small Tc However, mean- field theory neglects 

the infiuence of fluctuations, in particular, those linked to 
many-body processes resulting from interference between 
superconducting (SC) and density-wave (DW) pairings. 
Such a mix of pairing mechanisms is present at every order 
of perturbation theory for quasi-lD interacting electron 
systems and is responsible for the dynamical generation 
of d-wave superconducting pairing from spin fluctuations 
[laiMlIMlIlg]. For compounds like (TMTSF)2X, where 
superconductivity is just found in the close proximity of 
a SDW instability, the interference is expected to play an 
important role. Fully taken into account by the renor- 
malization group method (RG), the interference proved 
to have a non trivial impact not only on the mechanism 
and nature of the superconductivity [T4l - [T6l[30| . but also 
on the property of the metallic state above Tc. It carries 
over well below the crossover temperature scale for the 
onset of transverse coherent single electron motion where 
marked deviations from an ordinary Fermi liquid can be 
found [24l[29l[3ll|32]- The repercussions of these fluctu- 
ation effects on the stability of superconductivity under 
field are until now unexplored. 

In this letter we tackle the problem of upper critical field 
in quasi- ID superconductors close to a spin density- wave 
instability and clarify the influence of quasi-lD fluctua- 
tions on superconductivity under magnetic field. For this 
purpose, we extend the renormalization group (A'^-chains 
RG) technique j^l] to incorporate a finite Zeeman split- 
ting, and study the infiuence of interference on d-wavc sin- 
glet superconductivity under magnetic field. It is shown 
that the Hc2 exceed not only Hp, but also the mean- field 
prediction of H^2^^ to a great extent. The predicted non 
universal character of the effect, as a function of interac- 
tions and nesting alterations of the Fermi surface, corre- 
lates the degree of the interference between SC and DW. 
This give the possibility to confirm the relevance of the 
present theory by new experiments under pressure and 
field. 

Theory. — We consider a quasi-lD model of a square 
array of weakly coupling metallic chains, taken away from 
half-filling for simplicity. The hopping integrals ti obey 
the following anisotropic sequence 3> ^6 2> tc along 
the a-, 5-, and c-axis. The magnetic field is oriented such 
that its impact is restricted to the Zeeman splitting alone. 
Such conditions can be realized in practice when the field 



is oriented in the ab plane, which easily suppresses the 
orbital motion along the least conducting c axis at small 
enough t,. |25j . The remaining relevant electron spectrum 
in the remaining ah plane is given by 

£:±(k) = VF{±ka - kp) + Ukb) + fiBcrH, (1) 

where fiB is the Bohr magneton (in the following, we 
use h =: fiBH and h = 1, = 1), and ^b{kb) = 
—2tb cos kb — 2t'i^ cos 2kb. Here the longitudinal part of the 
spectrum has been linearized around ±kp, with the longi- 
tudinal Fermi velocity vp = 2taSiiikF. The next-nearest 
interchain hopping, tj,, breaks the perfect nesting condi- 
tion of the Fermi surface (FS). 

Electron-electron scattering matrix elements are ex- 
pressed in terms of interactions between electrons on oppo- 
site Fermi sheets |33]; namely as g\\ and gi± for backward 
scattering with parallel and antiparallel spins, and g2± 
for forward scattering. For the Hubbard model, we have 
9\\ =0, = ,92± = U/nvp = U, as bare normalized 
couplings, where U is the on-site repulsion of the Hub- 
bard model. We employ a momentum shell RG approach 
that is essentially equivalent to the one-particle irreducible 
scheme at one- loop level [T6l[24l[34l[35] . 

The one-loop flow equations so obtained are 

9,311 =311511 (4 -Ip)- 91^31^4", (2) 

digi± = - i9i±g2± + .92±5i±) {Ic + Ic) /2 

-(.9i±.9||+9||.9i±)(/p+/p")/2, (3) 

dmi^ = ,92±.92± (/p - Ic) - .9l±,9l±/c^ (4) 

where the momentum dependences of the couplings on the 
left-hand sides of Eqs. (E])-© are gi{kib, k2b, fcsb) with in- 
coming ki, — k2 and outgoing ka and — k4 = ki — k2 — ks. 
On the right-hand side, the momentum dependence of cou- 
pling products are in order: g(kib,kb — qi,,ki,)g{kf,,ki, — 
Qbiksb) for Cooper channel (particle-particle channel: SC 
fluctuations) that involves lQ^{kb,qb = ku — k2b)] and 
g{kib, Qb - kb, kb)g{kb, Qb - ksb, ksb) for the Peierls chan- 
nel (particle-hole channel: DW fluctuations) that involves 
Ip^{kb, Qb = k2b + kj,b)- The summations with respect to 
kb are taken as X^fe^i where N is the number of chains 
in the ah plane. The loop shell integration yields 

I^h^E^^ Qmi2 + \A^,'^\-E,/2) 

X <^ tanh — + tanh — — } , (5) 

\ AT 2T J ' ^ ^ 

where 

Ac'' = ^b{kb) - 6(9b - kb) - 6(fci6) + ^b{k2b) - ah, (6) 
A'^'^ = Cb{kb) + ib{kb - Qb) - Uk2b) ~ Uhb) - ah, 

(7) 

and 8(a;) = 1, 1/2 and for a; > 0, a; = and a: < 0, 
respectively. The scaled bandwidth is = Eqc^^ at the 
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RG step ^ > 0. The initial band width Eo{= 2wi?fci?) is 
taken as unity. The magnetic field acts as a cutoff for cer- 
tain logarithmic singularities through A"'' in the Cooper 
and Peierls channels, in a way similar to nesting deviations 
t'j^ for the Peierls singularity. 

From Refs. [231|33], the dimensionless response function 
for singlet type pairing at q = 0, reads 

9 /"OO 

xl^iT) = —j^M {zt;,%{i)Y (8) 
where the pair vertex obeys the flow equation 

diztCn (^s^n + 92±)/g")°'°. (9) 

The momentum dependences for the couplings and loop 
are {h, kf, - qi,, k'^) and /^'(gb, fc^)- denotes the 

even- (odd)-n-th Fourier component as (/(fcb, °^ = 
(/(fcf,, fc^) cos(sin)nA:b cos(sin)nA:^). In particular, (••■)§ is 
for s-wave and {•■•)% is for d-wave singlet, and (• • • )o is 
for p-wave and (• • • )i is for f-wave triplet pairings |14j . 

For the FFLO type pairing, q 0, and the formulation 
is modified due to the mixture of singlet and triplet pairing 
as follows: 

1 r<x> 

+ {4:^mficm, (10) 

where the pair vertex is renormalized as 

a _tie,o "'"/,tie,o^ rO I „4.te,o„ r4?i\e,o /I 1 N 

C^^ssg.k = - l^K^sq.n 92^-1 C +4sg,«5lJ-Jc In ■ (H) 

Similar definitions hold for the other pair response func- 
tions [241 . For the determination of the strongest singular- 
ity corresponding to the dominant instability, we compare 
the derivatives of the response functions, diXfj.- In the fol- 
lowing, we set tb = 0.1i?o, ij, = 0.01i?o, g\\ = 0.0 and finite 
repulsive U. 

Results. — Response functions. The magnetic field 
dependences of the derivative of the response functions, 
diXfi{H), for = 31 chains are shown in Fig. [TJ Here 
we set gi±_ = g2i_ = 1.0 {U = 1.0), and T = 10~^Eq. The 
SDW correlations are reduced (compared to d-wave) by 
nesting deviations due to t'j^. At low fields, the most dom- 
inant instability is for d-wave pairing with q = 0. (Here- 
after we call the non-FFLO type {q = 0) d-wave paring 
just as "d-wave pairing".) It is reduced for h > 10~'^i?o = 
r, reflecting the cut-off effect of h. The instability that 
closely follows is towards d-wave FFLO (dFFLO) order, 
which is also reduced for h > T, albeit less than d-wave, 
so that it dominates at high field. In the limit of very high 
field, h > IQ-'^Eo, the dFFLO and the transverse SDW 
(SDWjjj^) are prevalent with increasing H, feeding one an- 
other through the coupling g2_L as a result of interference. 




Fig. 1: Response functions, diXii{h), as a function of magnetic 
field for U = 1.0, T = 10"^£o and N = 31. The subscript of p 
and / denotes the Sz of triplet pairing; Sz = 1 is the equal spin 
pairing and = 1 is the triplet pairing with antiparallel spins. 
SDWz and SDW xy are the longitudinal- and the transverse 
spin density wave, respectively. 

H-T phase diagram. Next, we show the results of the 
temperature dependence of Hc2 for U =0.6, 0.8, and 1.0 in 
Fig. [2] Significant features are obtained. First the trace 
of Hc2{T), which is determined by the singularity in diXfi-, 
greatly exceeds not only the Pauli limit. Hp ~ 1.25Tc, but 
also the paramagnetic limit of the FFLO state, H^2^'^ in 
2D (~ l.TSTc) and 3D (~ 1.34rc ) isotropic cases 
Significantly enough, the limiting field ratio Hg^^'^ /T^ 
found here is non universal and depends on the strength 
of interactions. 

The origin of such an increase is primarily due to quasi- 
ID fluctuations introduced by SC-DW interference. The 
mixing between the Peierls and Cooper scattering channels 
in Eq. ^ induce an attractive coupling for d-wave super- 
conductivity from spin fluctuations, which is scale depen- 
dent all the way down to the lowest energy scale of the 
flow, that is Tc. This contrasts with the standard single 
channel or mean-field approach to superconductivity for 
which the attractive coupling that feeds the Cooper pair- 
ing channel is fixed and tied to the high-energy cutoff of 
the theory. This can be further checked by comparing the 
RG results with those of the RPA — ladder diagrammatic 
summation — limit, which is obtained by switching off the 
Peierls contribution to the flow of couplings constants in 
Eqs. (2)-(4) {I'^^=Q). There will be then no interference 
and by assuming an attractive coupling U < 0, an s-wave 
SC instability will thus be realized in the mean-field way, 
while keeping the same geometric condition on the FS. 

The results of ladder summation are also shown in 
Fig. [2] (the bottom solid curve). In this case Hc2{T) 
almost agrees with the ordinary one [3B|, Hc2{T) ~ 
Hpy^l — T/Tc around Tc, and extrapolates also towards 
Hp at zero temperature. As for the FFLO state in the 
ladder limit, H^.^^'~^{T) is higher than the isotropic value 
Hciau- This slight increase of H^^'^^iT) is due to a geo- 
metric property of the quasi-lD FS. In general, there is a 
geometric restriction to finite momentum FFLO pairing. 
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Fig. 2: 1^02 vs. temperature for U = 1.2, 1.0,0.8 (from top to 
bottom) for A*' = 15. The bottom curve is the Hc2 for U — 
—0.2 without the Peierls channels, i.e., no SC-DW interference 
(ladder diagrammatic summation). 



as shown in Fig. [3] (a). This restriction is partly removed 
when the FS shows nesting (Fig. [3] (b)). In the case of a 
quasi- ID FS, however, there is no restriction and FFLO 
pairing is enhanced, as shown in Fig. [3] (c). In the lad- 
der approximation, one finds that the traces for Hc2 /Tc vs 
T/Tc, as obtained for different interactions U (not shown), 
all scaled into the same universal curve at variance with 
the nonuniversality displayed by the full one-loop RG re- 
sults. 

Discussion. — Let us now discuss the implications 
of the present results for (TMTSF)2X compounds. For 
(TMTSF)2C104, the specific heat measurements of bulk 
critical field, H^2^^, along both a and b' directions agree 
with the Pauli limit at low temperature, but it is signif- 
icantly lower than H°2^'^^, as extracted from resistivity 
onset. This allows to define a rather large "transient re- 



gion", iJ, 



bulk 



, in the H-T plane where in- 



complete superconductivity is found. While there are no 
specific heat data available for (TMTSF)2PF6 under field, 
the values of H°2^^^ as obtained close to the critical pres- 
sure are similar to those of (TMTSF)2C104 at ambient 
pressure. This suggests that the transient region for both 
compounds is of similar size. 

The transient region observed experimentally should 
corresponds to the dFFLO region of the above two- 
dimensional theory (Fig. [5]). It seems to even exceed in 
size the one found experimentally at intermediate initial 
coupling U. However, dFFLO superconductivity is ex- 
pected to be quite sensitive to the finite impurity scatter- 
ing in these materials, an effect known to more pronounced 
for iJ;f I'LO [27] than for d-wave critical T^H = 0) [5H], 
which would result in some reduction of the transient re- 
gion. Furthermore, the influence of orbital pair breaking 
effect, not included in the present Zeeman field theory but 
which is finite in practice, would also contribute to the re- 
duction. 

Finally, a crucial test of the present theory of dFFLO 
state in materials like the Bechgaard salts would be clearly 
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Fig. 3: Possible FFLO pairings (q 7^ 0) for (a) isotropic, 
(b) 2D nested, and (c) quasi- ID FS. The FS at zero field is 
represented by the dashed lines. The pairs of (k, — k-|-q) locate 
on the FS split by the Zeeman field, (a) For the isotropic case, 
the pair of (k', — k' + q) cannot locate on the FS; the formation 
of the FFLO pair is difficult, (b) When the FS is nested, some 
of pairs like (k', — k' + q) can locate on the FS, but others like 
(k", — k"-|-q) cannot, (c) In the case of the quasi- ID FS, almost 
all pairs can locate on the FS, so that there is no geometrical 
restriction. 



the comparison of H-T phase diagram at various pres- 
sures, as it can be measured for instance by resistivity. 
In the present model the spin fluctuations induced attrac- 
tive coupling for d-wave superconductivity is known to be 
strongly pressure dependent via the anti-nesting param- 
eter t'f^, which besides U simulates pressure effects in the 
model [inillS]- Actually, the ratio Hc2/Tc decreases de- 
creases with increasing i'^ as shown in Fig. namely, the 
size of the transient region will shrink in size as pressure in- 
crease. The pressure dependent, non universal, reduction 
of Hc2 is a direct consequence of non universality linked 
to the interference between SDW and Cooper pairing. 

Summary. — In the present work, Hc2 for the quasi- 
ID superconductivity have been calculated for the Hub- 
bard model by the newly developed iV-chain RG under 
magnetic field. We have shown that when constructive 
interference between density-wave and Cooper pairings is 
fully taken into account, the Hc2 of a quasi-lD d-wave 
superconductor greatly exceeds not only the Pauli limit, 
but also the conventional paramagnetic limit of the FFLO 
state. In a material like (TMTSF)2C104, the results pre- 
sented here propose a new avenue for the reconciliation of 
transport [7J[5], magneto-torque [5], specific heat [TO] and 
NMR data [S]. The anomalous Hc2 enhancement found 
from this mechanism proves to be non universal in both 
interaction and anti-nesting strengths, allowing the possi- 
bility for a direct experimental test for the predicted cor- 
relation between the 7Jc2 -excess and Tc under pressure. 

In conclusion, we have proposed a mechanism of the 
anomalous broadening of the upper critical field region as 
seen from transport measurements in quasi-lD organic su- 
perconductors like the Bechgaard salts series. The mech- 
anism proposed is a direct consequence of the interplay 
between magnetism and unconventional superconductiv- 
ity in these materials under magnetic field. 
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Fig. 4: Upper critical field for different anti-nesting paramters 
[t'i/Eo = 0.005,0.01,0.02) for tt/Eo = 0.1, U = 0.8 and iV = 
15. The effect of t'^, mimics the pressure effect in (TMTSF)2X. 
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